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We present a study of the fiber bundle model using equal load sharing dynamics where the breaking 
thresholds of the fibers are drawn randomly from a power law distribution of the form p{b) ~ b _1 
in the range 10 _/3 to 10 /3 . Tuning the value of f) continuously over a wide range, the critical 
behavior of the fiber bundle has been studied both analytically as well as numerically. Our results 
are: (i) The critical load o c ([3,N) for the bundle of size N approaches its asymptotic value cr c (/i) as 
a c (/3,N) = o c (/3)+AN~ 1 / I/ (^ where <r c (/3) has been obtained analytically as cr c (/3) = 10'®/(2/3e In 10) 
for /3 > fl u = 1/(2 In 10), and for (3 < j3 u the weakest fiber failure leads to the catastrophic breakdown 
of the entire fiber bundle, similar to brittle materials, leading to cr c (/3) = 10~^; (ii) the fraction of 
broken fibers right before the complete breakdown of the bundle has the form 1 — 1 / (2/3 In 10); (iii) 
the distribution D( A) of the avalanches of size A follows a power law D( A) ~ A”' with £ = 5/2 for 
A A c (/3) and f = 3/2 for A -C A c (/9), where the crossover avalanche size A c (/3) = 2/(1—elCP 2 ^) 2 . 

PACS numbers: 64.60.Ht 62.20.M- 02.50.-r 05.40.-a 


1. INTRODUCTION 


Natural disasters like land slide, mine collapse, earth¬ 
quake cause great losses in human lives and property. It is 
therefore primarily important to understand the underly¬ 
ing mechanisms of the failure processes so that the losses 
can be minimized by providing a precursor. Similarly for 
engineers the strength of material is a major quantity in 
order to make huge constructions like bridges, buildings 
etc. Due to these standing requirements, during the last 
two decades, huge amounts of scientific efforts have been 
invested to explore the microscopic mechanism and rup¬ 
ture process of disordered materials. It has been revealed 
that the disorder plays a crucial role in determining the 
strength of material and also in the fracturing process 

ffi. 

Models of materials in the form of a bundle consist¬ 
ing of a large number of parallel massless elastic fibers 
are well known to be simple examples of critical sys¬ 
tems exhibiting non-trivial breakdown properties [lj-[5j. 
These systems are called the Fiber Bundle Models (FBM) 
where individual fibers have randomly distributed break¬ 
ing thresholds. Typically, on increasing the externally 
applied load a per fiber, the entire fiber bundle fails at a 
critical load a c per fiber. It is also known that for a < a c , 
larger the external load, more extensive is the response of 
the system in terms of the number of fiber failures. This 
number diverges as a —► a c from below and for er beyond 
a c all fibers eventually fail with certainty. Therefore, a c 
is looked upon as the transition point from a local to the 
global failure of the bundle @]. 

In the fiber bundle model a set of N parallel fibers is 
clamped at one end and an external load is applied at 
the other end 3, [sj. Every fiber i has its own break¬ 
ing threshold 6,. If the tensile stress acting through it 
exceeds bi it breaks. Random numbers {6^} are drawn 
from a probability distribution p(b) and they are assigned 
as the breaking thresholds of the individual fibers whose 


cumulative distribution is P(b) = f 0 p(z)dz. 

In many FBMs, stress is treated as a conserved quan¬ 
tity. During the failure of an individual fiber the stress 
is released and it gets distributed among the remaining 
intact fibers. Depending on how the released stress is 
distributed among the intact fibers there exists various 
models in the literature. Among these FBMs the Equal 
Load Sharing (ELS) model is the most well known [8j-|l0]. 
Here the released stress is distributed equally among all 
the remaining intact fibers. Most of the results of this 
model have been calculated analytically and also this 
model is computationally easier to tackle. On the other 
hand in the Local Load Sharing (LLS) model, the re¬ 
leased stress is distributed equally only to the nearest 
surviving neighbors mua. In the LLS model, most 
of the results have been obtained numerically. A fiber is 
strained when some amount of stress acts through it. For 
ELS, the clamps at two ends of the bundle may be treated 
as infinitely stiff and therefore under a certain applied 
load, all fibers are strained by equal amounts and con¬ 
sequently the magnitudes of the stresses acting through 
the intact fibers are also equal. On the other hand, if 
the clamps are elastic, different fibers are strained differ¬ 
ently and their stress values are also different as is the 
situation in the LLS model. A third model, intermedi¬ 
ate between ELS and LLS, has also been considered in 
the following way. In this model the released stress is dis¬ 
tributed non-uniformly and the share amount received by 
an intact fiber depends inversely to some power of the dis¬ 
tance of separation from the broken fiber fbt . A number 
of other processes have been studied in the framework 
of fiber bundle models. For example, how the damage 
evolves due to an environmentally assisted a ging process 
in a fiber bundle model has been studied in 14]. In this 
paper, we study the breakdown properties of the fiber 
bundles with ELS dynamics. 

Let a be the uniform applied load per fiber initially 
when all fibers are intact. The total amount of external 



2 


load is then F = Na. This externally applied load gets 
distributed within the bundle in a series of T successive 
time steps. Let us denote x t as the stress per intact fiber 
after t-tii relaxation step. Since more and more fibers 
break, the stresses acting through the remaining intact 
fibers increase. When a is the applied load per fiber, all 
fibers with bi < a break. This stress is now distributed to 
N[1 — P{ct)} intact fibers on the average. After the first 
step if x\ is the stress per fiber then F = 7Vxi[l — P(cr)]. 
Consideration of the same mechanism in successive steps 
one can write: 


F = N Xl [l-P{a)] = Nx 2 [l-P{x 1 )} = Nx 3 [l-P{x 2 )}.. 

( 1 ) 

This process terminates after T steps when the amount 
of stress released is not sufficient to create further failure 
of fibers. 

If x is the applied load per intact fiber in the stable 
state, then one can write the external load F(x) as a 
function of x. which is F(x) = Nx[l — P(x)] @,[l5|. For a 
specific value of x = x c . F(x) is maximum which suggests 
the condition: 1 — P(x c ) — x c p(x c ) = 0. For example, for 
a uniform distribution of breaking thresholds one gets 

<r c = F c /N =1/4 0. 


The failure properties of materials are highly depen¬ 
dent on the extent of disorder inherent in them. In the 
FBMs, this disorder appears in the breaking thresholds 
of the individual fibers. In this regard, the power law dis¬ 
tribution of breaking thresholds is an extreme case of het¬ 
erogeneous disorder, where a large number of fibers have 
very small breaking thresholds, their numbers decreases 
as breaking thresholds are increased, leading to few fibers 
with large breaking thresholds. It is already known in the 
literature that the probability of getting a warning of im¬ 
minent breakdown of the system is higher when the ma¬ 
terial is more heterogeneous [16] . Such cases of extremely 
heterogeneous disorder has not been very well studied in 
the literature of fiber bundle model. Another form of 
strong heterogeneity has been studied where a fraction 
of fibers are completely unbreakable and the breaking 
thresholds of the rest are drawn from some distribution 


13,118]. In this paper we therefore address the problem of 


FBMs with highly heterogeneous power law distributed 
breaking thresholds of individual fibers. As a first step 
we study the simpler problem of ELS dynamics in this 
model, study of the LLS version will be taken up in a 
future publication. 

FBMs have a wide variety of applications. It is a ver¬ 
satile tool to understand conceptually the underlying mi¬ 
croscopic mechanism of fatigue II , failure of composite 
materials {20|, landslides [21] etc. Moreover, the ELS 
version of the FBM studied here may be used to study 
the traffic jams in roads j22j|. The traffic flow capacities 
of the roads can be mapped to the breaking thresholds 
of individual fibers. The highly disorder flow rates may 
occur in a traffic network with few highways and a large 



FIG. 1: The initial external load <t c (/ 3) per fiber given in Eqn. 
5 (solid line) matches excellently with its numerical estimates 
(open circles). 


number of narrow roads connecting the highways. 

In section 2 we describe our study of the fiber bundle 
model with power law distributed breaking thresholds. 
We also describe different analytically obtained results 
characterizing this bundle and their numerical supports. 
In section 3 the statistics of avalanche size distribution 
have been described. We summarize in section 4. 


2. HIGHLY DISORDERED FIBER BUNDLES 


In this paper we report the results of our study of the 
breakdown properties of a fiber bundle where the break¬ 
ing thresholds of the individual fibers are power law dis¬ 
tributed. As in other FBMs, the only source of disor¬ 
der in our model is the random distribution of breaking 
thresholds. Therefore, the individual breaking thresholds 
bi are drawn from a probability distribution p(b) ~ 6 _T 
with 7 = 1 . Initially N uniformly distributed random 
numbers qt are drawn within — 1 < qi < 1 and the break¬ 
ing threshold bi = 10 ^ 9i for the i-th fiber is assigned. 
Consequently, the probability distribution takes the form 
p(b) ~ 6 _1 within the range 10 _/3 to 10 ^ 161. 

Here, we use the same formulation described in sec¬ 
tion 1 to obtain the breaking strength of the bundle 
cr c as a function of the cut-off parameter /3 when the 
breaking thresholds {bi} are power law distributed. The 
constant of proportionality can be evaluated from the 
normalization condition, which gives the functional form 
p(b) = 6 - 1 /(2/llnl0). As a result, the cumulative prob¬ 
ability distribution is given by, 

P(b) = [ p(z)dz = lni>/(2/31nl0) + 1/2. (2) 

Jw-P 

In this case we obtain the expression of F{x) as 


F{x) = Nx[l/2 — lnx/(2/31nl0)]. 


(3) 


Clearly the function F(x) has a maximum at x = x c 
for which dF(x)/dx = 0. This yields x c = 10^/e 
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FIG. 2: Variation of the critical load a c (/3,N) on the system 
size TV for /? = 0.225 has been exhibited. Plot of a c (/3,N) — 
<Tc{P) vs. TV -0 ' 624 with a c (/3) = 0.596 shows a nice straight 
line that passes very close to the origin. 


and the total critical applied load is F c = F(x c ) = 
TV10^/(2/3e In 10). Thus the critical initial applied load 
per fiber is given by 

a c {/3) = F c /N = 10' 9 /(2/3e In 10). (4) 

Let b* denote the minimum of the breaking thresholds. 
Since the definition of x c signifies that a bundle fails com¬ 
pletely at this point then, the condition b* = x c i.e., 
10 _/3 = 10 9 /e fixes the upper bound of p denoted as 
Pu = 1/(2 In 10 ) for which the weakest fiber failure leads 
to the complete breakdown of the bundle. Thus we have 
the complete expression for a c (/3): 




FIG. 3: (Color online) (a) Plot of v(j3, TV) vs /3 for systems of 
different sizes. The value of v(/3, N) calculated using the four 
bundle sizes from N = 2 10 to 2 16 (black), 2 14 to 2 20 (blue) 
and 2 18 to 2 24 (red); TV is increased from left to right, (b) 
A collapse of the data of the same three system sizes works 
excellent when the P axis has been suitably scaled. 


10 ^/( 2 e In 10 9 ) for P > j3 u 

a (5) 

10 1 for P < p u 

The above expression for the average critical applied load 
per fiber <j c (P) for a given value of cut-off parameter P is 
valid only for infinitely large bundles, i.e., TV —> oo. 

The width of the distribution of breaking thresholds 
increases with p and the critical threshold <r c (P) varies 
accordingly. For P = 0, all fibers have the same break¬ 
ing thresholds equal to unity and therefore cr c ( 0 ) = 1 . 
When P is small, the minimum breaking threshold is 
high enough, and very close to unity. When the external 
stress per fiber is raised to reach the minimum breaking 
threshold, it breaks. The released stress is distributed 
among the remaining fibers and is sufficient to break all 
other fibers. This mechanism, when failure of the weak¬ 
est fiber ensures the global failure of the entire bundle 
is analogous to the brittle fracture. This situation con¬ 
tinues till p reaches p u and therefore a c {P) decreases as 
the strength of the weakest fiber, i.e., 10 - ' 3 . When p in¬ 
creases further, gradually fibers of high breaking thresh¬ 
olds appear and they take over the control. Consequently, 
( t c (P ) must increase with P for large p with a minimum 
at p = p m . The value of p m is obtained using the condi¬ 
tion da c {p)/dp = 0 in Eqn. (5) at p m = l/(ln 10), which 
is twice the value of p u . 


Numerically cr c (p) is obtained in the following way. 
For a given value of p we first calculate the critical load 
per fiber er“(/?,TV) for a particular fiber bundle a having 
TV fibers with a given set of breaking thresholds {bi}. 
This calculation is repeated over a large number of un¬ 
correlated bundles a and their critical loads are averaged 
to obtain a c (P,N) = (cr“(/3, TV)}. The entire calculation 
is then repeated for different values of TV. 

To obtain <r“(/3, TV), the breaking thresholds are ar¬ 
ranged in increasing order ( b < bfy < b < ... < 
b? N ^). The bundle will support the initially applied load 
per fiber (a) if a < b or ctTV/(TV — 1) < b “ 2 ^ or 
trTV/(TV — 2) < or ... crTV < b? N y If all these in¬ 
equalities fail to satisfy then the bundle will no longer 
support the load, it will break apart. Now if a is such 
that it is sufficient to break n fibers, then at this stage 
the bundle will support the load if trTV/(TV — n) < fr“ n+1 ) 
i.e., 

a<[(N-n)/N]b? n+iy (6) 

The term in the parenthesis of Eqn. ( 6 ) decreases with 
n and 6“ n+1 j is an increasing function of n as thresh¬ 
olds are arranged in increasing order. So, the function 
at the right hand side of Eqn. ( 6 ) has a maximum at 
some n and if the a is raised at this maximum value, 





4 



FIG. 4: Plot of the fraction of broken fibers fb(fi) for a par¬ 
ticular value of p right before the complete breakdown of the 
bundle is plotted for the analytical expression given in Eqn. 
9 with solid line. The numerically obtained data represented 
by open circles matches very well with the analytical curve. 

the bundle will break immediately. So the maximum of 
[(TV — n)/AT]&“ n+1 j determines the critical load per fiber 
for the bundle a. Therefore HH, 



io° to 1 nr a io 3 io J io s io 6 



°ciP> N ) = max {^(lp 


TV -1 TV - 2 1 

h a h a h a 

~N~ b ( 2 )' ~N~ b W’ N b W. 


(7) 

We now assume that the average value of the critical 
load per fiber cr c (/3, TV) for a given value of j3 and for the 
bundle of size TV converges to a specific value <r c (/T) as 
TV —> oo according to the following form: 


a c (/3,N)-a c (P) = AN~ 1 ^ (8) 


FIG. 5: (Color online) (a) Log-log plot of the binned data 
for avalanche size distribution D( A) vs A for /? = p u = 
1/(2In 10) for TV = 2 16 , TV = 2 18 ... 2 26 (from left to right), 
(b) A finite-size scaling works well: D(A)N rl against ATV” ? 
exhibits a good collapse of data with 77 = 1.007 and £ = 0.671 
implying £ = 77 /C = 1-50(1). This value is consistent with 
the directly measured value of 1.50(2) from the slopes in the 
intermediate region. The crossover is not observed here since 
A c = 00 for this particular value of /3. 


where i/(/3) is a critical exponent for the cut-off parameter 
[3. We have plotted cx c (/3, TV) against TV’ _1 /*'0 3 ) for TV = 
2 18 to 2 24 , TV being increased by a factor of 4 at each 
stage. For a particular value of /3 we have used different 
trial values of v{fi) so that for a specific value of is(/3) 
the plot fits (by least square fit) to the best straight line. 
Using this best value of i/(/3) and on extrapolation to 
TV —> 00 we obtained er c (/3). In Fig. 1 we have exhibited 
an excellent matching of the analytical and the numerical 
values of er c (/3) for the range 0 < /3 < 2. 

We now investigate the dependence of the finite size 
correction exponent i/(/3) on the cut-off parameter /3. We 
recall that in the case of a uniform breaking threshold 
distribution, the plot of <t c (TV) — cr c as a function of TV” 1 / 1 ' 
gives an excellent straight line with a c = 1/4 and v = 3/2 


24-27]. Similarly, for our model of highly disordered 


FBM, the plot of cr c (/3,TV) — cr c (/3) against is 

carried out for different values of ft. For example, we 
obtain the best possible value of v(j3) to be 1.603 for f} = 
0.225 shown in Fig. 2. In this way the critical exponent 
z/(/3) is calculated for different f3 and its variation is shown 
in Fig. 3(a) using TV = 2 10 to 2 16 , 2 14 to 2 20 and 2 18 to 
2 24 . The value of iz(/3) first increases, attains a maximum 
value ss 1.63, then decreases and saturates to 1.5 with 


further increment of /?. The same data in Fig. 3(a) when 
plotted against (/3 — /3 U )TV°- 33 shows a good collapse as 
shown in Fig. 3(b). Thus we conclude that the curve for 
v{fi) retains its nature for large bundle sizes. 

Next we calculate the fraction of broken fibers fb{/3) 
just before complete breakdown of the bundle as a func¬ 
tion of the cut-off parameter fj. Since at x c the fiber 
bundle fails completely, so the quantity /&(/3) is calcu¬ 
lated as: 

fb(/3) = f p(x)dx = 1 - 1/(2^ In 10). (9) 

Jio-e 

As the fraction of broken fibers fb{P) is a positive quan¬ 
tity thus the condition 1 — l/(2/31nl0) > 0 again re¬ 
produces the result that for f) < 1/(2 In 10) the weakest 
element failure leads to the catastrophic breakdown of 
the bundle. 

Numerically /(,(/?), for a given value of /3 is calculated 
in the same way as described previously in the case of 
cr c {/3). The external load is increased quasi-statically un¬ 
til the bundle fails. Just before complete breakdown of 
the bundle, the fraction of broken fibers is calculated for 
a particular TV and averaged over large number of sam- 
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where, 
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10 10 A 10 10 10 ' 10 


r(x) = 1 — 


xp(x) 

1 - P(x)' 


( 12 ) 


The expression for D(A) can be simplified to the follow¬ 
ing form |29}: 


£>(A) _ A a 2 e A p{x c ) —a/Ac^ 

N A! |r'0r c )| l ' 


(13) 


with 


FIG. 6: (Color online) The avalanche size distribution for f3 
= 0.22 (black), 0.24 (red) and 0.28 (blue) (from right to left) 
for bundles of size N = 2 24 . Slopes of the curve are ss 1.5 and 
Ri 2.5 for small and large avalanche sizes. The crossover size 
A c (/3) are approximately 11741, 200.4 and 31.66 respectively 
evaluated using Eqn. 17. 

pies to obtain fb(/3,N). Then this procedure is repeated 
for six values of N = 2 , 2 18 ,..., 2 26 and an extrapolation 
on fb(/3, N) as N —> 00 yields /&(/3). In Fig. 4 the numer¬ 
ically obtained results are compared with the analytical 
one, indicating a good agreement. 


c r'(x c ) 2 (x c — b*) 2 ' 1 j 

Using the Stirling approximation A! = A A e _A V27rA, 
Eqn. (13) can be written as 

= CA- 5 / 2 (l-e- A/Ac ), (15) 

Where C = {2i:)~ 1 ^ 2 p{x c )/\r'(x c )\ is a constant. From 
Eqn. (15), a clear evidence of crossover in the exponent £ 
around the avalanche size A c is prominent. So we have: 


3. AVALANCHE SIZE DISTRIBUTION 


D( A) f A 3/2 for A < A c , 

-OC \ 

N [ A -5 / 2 for A > A c . 


(16) 


In a stable fiber bundle, the stress acting through ev¬ 
ery intact fiber is less than its breaking threshold. Now, 
if the externally applied load is suitably raised so that it 
becomes equal to the breaking threshold of the weakest 
fiber, then this fiber breaks. This triggers a cascade of 
fiber failures which finally ends when the bundle attains 
a new stable state. The total number A of fibers that 
fail in this event is called the avalanche size. Starting 
from a completely intact fiber bundle the global failure 
of the entire bundle may be attained by raising the ex¬ 
ternal load in such a quasi-static process, causing a se¬ 
ries of avalanches. The probability distribution D( A) 
is regarded as an interesting quantity to study. It is 
well known that for the uniform distribution of breaking 
thresholds of individual fibers and for the ELS dynamics, 
the probability distribution is a power law [28] 

D( A)~A"«, (10) 

with £ = 5/2. In the following we would see that in 
our case of power law distributed breaking thresholds the 
exponent £ undergoes a crossover from 3/2 to 5/2. 

To exhibit the crossover behavior we follow the method 
in [ 29 }. For a bundle having large number of fibers, the 
number of avalanches of size A is given by [28( | 

= A A , 6 - p{x)r(x )[1 - r(x)] A ~ 1 e Ar ^dx, 

( 11 ) 


In our case, we use power law distribution p(b) ~ 5 _1 
in the range from lO”' 3 to 10' 3 to obtain r'(x c ) = —e/10' 3 , 
x c = 10 13 /e and b* = 10~' 3 . Substituting these values in 
Eqn. (14) we get the crossover avalanche size: 

= (1-elO- 2 / 3 ) 2 ' (17) 

This crossover phenomenon has also been studied using 
numerical simulations. For (3 = 1/(2 In 10) Eqn. (17) 
yields A c = 00 . Thus only the £ = 3/2 power law is 
observed as any avalanche of finite size A is less than the 
value of A c at this particular value of /3. In Fig. 5(a), 
the numerical data for the avalanche size distribution for 
/? = 1/(2 In 10) has been plotted for six different values of 
N starting from N = 2 16 to 2 26 ; N being increased by a 
factor of 4 at each stage. For N = 2 16 to 2 22 the data has 
been averaged over 10 6 samples and 400000 and 100000 
samples for 2 24 and 2 26 respectively. A finite-size scaling 
has also been done in Fig. 5(b) by use of suitable powers 
of the bundle size N. This indeed exhibits an excellent 
data collapse confirming the following scaling form: 

D(A)N r> ~ g[A/N c ] (18) 

where Q(y) is an universal scaling function of the scaled 
variable y = A/N^. The best possible tuned values of the 
scaling exponents obtained are r) = 1.007 and £ = 0.671. 
Using these scaling exponents the value of £ = rj/( = 
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FIG. 7: Plot of the average number of avalanches A(IV) re¬ 
quired to break the bundle of size N on a log-log scale: 
for /3 = 1/(2 In 10), A{N) ~ A 0 ' 337 (filled circles) and for 
P = 0.240, A (N) ~ AT 0 - 985 (open circles). 

1.50(1) is calculated, which is a very well tally with the 
analytical result of 3/2. 

We have also tried the same analysis for p = 0.22,0.24 
and 0.28. Using Eqn. (17) we have obtained A c (/3) = 
11741,200.4 and 31.66 respectively. A clear evidence of 
the crossover in the exponent £ around A = A c (/3) is 
observed as shown in Fig. 6 for N = 2 24 . The slope 
of the curve gradually crosses over from ss 1.5 to ~ 2.5 
for large values of A. It has also been observed that as 
P is increased, A c (/3) gradually shifts towards the origin 
and therefore the regime over which £ = 5/2 is valid, 
gets extended. Such a crossover has been observed ear¬ 
lier in |29j, [30] for the FBM with uniform distribution of 
breaking thresholds ranged between a certain lower cut¬ 
off bi c and unity. Here, avalanche sizes smaller (larger) 
than some crossover size A (bi c ) correspond to avalanche 
size exponents 3/2 (5/2). This implies that in our model, 
even for the highly heterogeneous distribution of break¬ 
ing thresholds, similar crossover between the same two 
exponents takes place across the crossover avalanche size 

A c 00). 

It has also been observed that the total number of 
avalanches A (N) depends on the system sizes N as N x , 
where x = 0.336 and 0.985 for /? = 1/(2In 10) and 0.240 
respectively. The log —log plot of A(N) against N for 
these two values of p fits to excellent straight lines as 
shown in Fig. 7. We conjecture that x nray be 1/3 and 
1 exactly for P = (3 U and P > p u respectively. 

4. SUMMARY 

Properties of the fiber bundle model have been stud¬ 
ied using equal load sharing dynamics where the breaking 
thresholds of the fibers have been assigned from a power 
law distribution p(b) ~ b~ 1 in the range from 1CV' 3 to 
10 /3 . Variations of different quantities characterizing the 
bundle have been studied with the cut-off exponent p. 


The critical external load per fiber a c (P) required for the 
global breakdown of the bundle as well as the fraction 
fb(P) of broken fibers right before it are estimated both 
analytically as well as numerically, and a good correspon¬ 
dence has been observed. For very small and very high 
values of p the breaking strength of only a single fiber 
determines the critical strength of the entire bundle. For 
example, for very small /3, it is enough to tune the ex¬ 
ternal load to the strength of the weakest fiber which 
then triggers a large avalanche and the entire bundle fails, 
implying <r c (/3) = 10 -/3 . Such a behavior continues till 
P = [3 U and this regime is analogous to the brittle fail¬ 
ure of materials. When P is raised beyond P u , equating 
the external load to the strength of the weakest fiber is 
no more sufficient for the global failure, the large num¬ 
ber of fibers with breaking thresholds near the minimum 
still dominate. Consequently the number of avalanches 
required for the breakdown of the bundle gradually in¬ 
creases and the <r c (/3) slowly increases from the weakest 
strength of IO - ' 3 but for P > p u , cr c (P) remains smaller 
than 1CV ,3 “. Therefore, a minimum in a c (P) is reached 
at p m = 2P u and from this point, a c (p) starts increas¬ 
ing. As a result, cr c {P) becomes equal to 10 _/3 " again 
and then increases indefinitely. For very large P the ex¬ 
ternal load must be raised to cr c (P) ~ 10^ to break the 
strongest fiber of breaking threshold around 10' 3 . This 
salient feature is a direct consequence of the power law 
distribution of the breaking thresholds. 

More interestingly, we have also shown numerically 
that the critical load a C (P, N ) approaches its asymptotic 
value as cr c {P,N) = a c (p) + AN The finite size 
correction exponent zz(/3) is first seen to increase sharply 
with P , reaches a maximum, then decreases and finally 
converges to a value ~ 3/2. Statistical analysis of the 
avalanche sizes have been done. The avalanche size dis¬ 
tribution follows a power law and the associated expo¬ 
nent £ crosses over from 3/2 to 5/2 through a crossover 
avalanche size A c (/3). 
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